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OO ' Precision measurement of the stop mass at the ILC is done in a method based on 

. cross-sections measurements at two different center-of-mass energies. This allows to 

minimize both the statistical and systematic errors. In the framework of the MSSM, 
a light stop, compatible with electro-weak baryogenesis, is studied in its decay into a 
charm jet and neutralino, the Lightest Supersymmetric Particle(LSP), as a candidate 
of dark matter. This takes place for a small stop-neutralino mass difference. 
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1 Introduction 



In this study wc aim at the minimisation of the systematic uncertainties and of the statistical 
' error [T]. This is achieved by using a method which allows to increase the precision in two 

OO . ways. We deal with a ratio of cross-sections at two energy points. This takes care of the 

CO ■ systematic uncertainties by cancelations and we choose one of the energies to be at the 

threshold were the sensitivity to mass is maximale. We will show that even though we are 
' dealing with more realistic data than in [2], we improve substantially the precision in the 

mass measurement. As in [2], we are considering the MSSM with R Parity conservation 
and a scenario in which a light stop co-annihilates with the Lightest Supersymetric Particle 
(LSP), the neutralino, to produce the right amount of dark matter relic density, namely, 
within the experimental precision of WMAP and the Sloan digital sky survey [3] . Together 
^ , with a light Higgs, a light right-handed stop also supports electroweak baryogenesis. Our 

^ ' data now include hadronization and fragmentation of the stop before its decay as well as 

5^ , fragmentation of the charm of the decay. This provides a rather big smearing of the particles 

produced and together with gluon radiation increases the number of jets. We will use two 
different approaches. First wc will optimize a set of sequential cuts as in [5], then we will be 
using a multi-variable optimization of the neural-network type, the Iterative Discrimination 
Analysis(IDA). We do take also advantage of the polarization since we deal with an almost 
right-handed stop as required for E.W. baryogenesis. This allows us to enhance the signal 
while getting rid of a big part of the background. 

2 Mass Precision Measurement: the Method 

• The production cross-section of stop pairs e^e^ — > ti 4 is represented to next to 
leading order (NLO), as a function of the energy for two hypothetical values of the 
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Figure 1: Precision in Pair Production Cross-Section 

stop mass, 122. 5 GeV and 123.5 GeV, shown in Figure [TJ 

• In the lower left figure the seale has been blown up and one can see that the sensitivity 
to small mass difference is high at or close to threshold while in the lower right figure 
one sees that it is not the case at peak value. 

• We will define a parameter Y, as a ratio of production cross-sections at two energy 
points. This will reduce the systematic uncertainties in Y from the efficiencies as well 
as from the beam luminosity measurements between the two energy points. 

• One of the energy points is chosen at or close to the production energy threshold. This 
provides an increased sensitivity of Y to mass changes. 
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Figure 2: Precision in Determination of the Stop Mass 

Y{Mx ~ = '^iVs^)<^thLth 

a is the cross-section in [fb] , N the number of detected data, B is the number of estimated 
background events, s is the square of the center of mass energy, e the total efficiency and 
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acceptance and L is the integrated luminosity. The suffix (th) is used for the point at energy 
threshold and (pk) for the energy peak. M^: is the mass to be determined with high precision. 

As an example, one assumes 3% precision for Y, The blue hashed region represents our 
measurements. One obtains a precision AM^, ±0.2%, the 2 vertical arrows. 

In the method, we determine the stop mass by comparing Y with the theoretical calcula- 
tion of the cross-sections to next to the leading order (NLO) for both QCD and QED.It has 
been done for +80% polarizations for the e~ beam and —60% polarization for the e+bcam. 

3 The Channel Studied e+e" ^ii~A^ cX^cX^ 

A scan in the super-symmetry parameter space [5] has shown that a stop mass of 122.5 
GeV and a neutralino mass of 107.2 GeV are consistent with baryogenesis and dark matter. 
The process and the background channels are listed below with their cross-sections with and 
without polarization. 



Process 


Cross- 


section [pb] at ^/s 


= 260 GeV 


Cross 


section [pb] at y's 


= 500 GeV 


P(e-)/P(e+) 


0/0 


-80%/+60% -1 


-80%/-60% 


0/0 


-80%/-h60% -t 


-80%/-60% 


tit\ 


0.032 


0.017 


0.077 


0.118 


0.072 


0.276 


w+w- 


16.9 


48.6 


1.77 


8.6 


24.5 


0.77 


zz 


1.12 


2.28 


0.99 


0.49 


1.02 


0.44 


Wev 


1.73 


3.04 


0.50 


6.14 


10.6 


1.82 


eeZ 


5.1 


6.0 


4.3 


7.5 


8.5 


6.2 


qq, q^t 


49.5 


92.7 


53.1 


13.1 


25.4 


14.9 


tt 


0.0 


0.0 


0.0 


0.55 


1.13 


0.50 


2-photoii 

Pt > 5 GeV 


786 






936 







Table 1: The Cross-sections at ^/s = 260 GeV and y/s = 500 GeV for the signal and 
Standard Model baekground are given for different polarization combinations. The signal 
is given for a stop mixing angle of 0.01 and for a stop of mi — 122.5 Gev,consistant with 
E.W. baryogenesis. The e~ negative polarization values refer to left-handed polarization 
and positive values to right-handed polarization. 



3.1 Simulations Characteristics 

The signal and background channels were generated with Pythia(6.129), the simulator 
Simdet(4.03) and for the beamstrahlung Circe(1.0)[6]. They were generated in proportion 
with their cross-sections. 

• Hadronization of the ti quark and the fragmentation of the charm quark come from 
the Lund string fragmentation model. We use Peterson fragmentation [7] . 

• The stop Hadronization and fragmentation are simulated using T. Sjostrands code as 
described in detail by A.C.KraanfJ. The stop quark is set stable until after fragmen- 
tation, then it is allowed to decay. The stop fragmentation parameter is set relative 
to the bottom fragmentation parameter = efcrn^/m? and et = —0.0050 ± 0.0015. 
Later improvements at LEP and a factor two improvement assumed at ILC leads 
to Aet=0.6xl0~6, as detailed in 0. The charm fragmentation is set from LEP to 
ec = -0.031 ±0.011. 



LCWS/ILC2007-June2,C.Milstene. 



• The mean jet multiplicity increased for the data with fragmentation included. 
4 The Analysis 

The ntuple analysis code [5] which incorporates the Durham jet algorithm is used. The pre- 
selection and selection cuts are discribed in detail at both energies in [8] . For a right-handed 
stop, a run with right-chirality is favorable as shown in Table [l] one expects at 500 GeV a 
luminosity of 200 fb^l out of the 500 fb^l. The luminosity of 200 fb^l is used for the 0/0 
beam polarization as well as a point of comparison. 

4.1 The sequential cuts 

Were made as similar as possible at the two energies to aim at the cancellation in Y of the 
systematics. The cuts and their detailed results are given in [5]. In this analysis we allow 
two, three or four jets with the request that Eni25 GeV for the lowest energy jets. The 
charm tagging is extracted using the ZVTOP software. The product of Charm tagging of 
the two jets with the biggest charm probability is used to separate the signal from its main 
background, since the Wez^ has at most one charm jet, whereas the signal has two charm- 
jets. The backgrounds and the signal efficiencies are shown at the two energies after the 
sequential cuts in Tabled] 

4.2 Iterative Discriminant Analysis (IDA) 

Combines the kinematic variables in parallel. The same kinematical variables and simulated 
events are used than in the cut-based analysis. A non-linear discriminant function followed 
by iterations enhances the separation signal-background. Both signal and background have 




Efficiency Efficiency 



Figure 3: Detection Efficiency and Background Events at 500 GeV(left) and 260 
GeV (right). 

been divided in equally sized samples, one used for the training, the other as data. We will 
make two IDA iterations in our final analysis |S]. The results are shown after a first IDA 
iteration for which one keeps 99.5% of the signal efficiency followed by a second iteration. 
We assume the same luminosities and polarizarions than for the sequential based analysis. 
The results are in Table O The events < N show the number of events corresponding to a 
single event. 
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Sequential Cuts 


IDA 


IDA 


Pe-/Pe+ 


generated 


0/0 


Pol+ 


generated 


0/0 


Pol+ 


0/0 


Pol+ 


0/0 


Pol+ 


tit* 


50,000 


544 


1309 


50,000 


5170 


12093 


619 


1489 


9815 


22958 


W+W~ 


180,000 


38 


4 


210,000 


16 


2 


11 


1 


<8 


< 1 


ZZ 


30,000 


8 


7 


30,000 


36 


32 


<2 


<2 


20 


18 


Weu 


210,000 


208 


60 


210,000 


7416 


2198 


68 


20 


1719 


510 


eeZ 


210,000 


2 


2 


210,000 


<7 


<6 


3 


2 


<7 


<6 




350,000 


42 


45 


350,000 


15 


17 


16 


17 


18 


21 


tt 








180,000 


7 


7 






1 


1 


2-pliotoii 


1.6x10^' 


53 


53 


8.5 X lO'' 


12 


12 


27 


27 


294 


294 


Total Background 


351 


171 




7509 


2274 


127 


69 


2067 


851 


S/B 


1.5 


7.6 




0.7 


5.3 


4.9 


22 


4.7 27 


Efficiency 


0.340 


0.212 


0.387 


0.416 



Table 2: Signal and background generated to NLO and after selection cuts are shown at 
y/s = 260 GcV and 500 GcV, for total luminosities of 50 fb~^ and 200 fb~^ , respectively 
and the signal efficiencies. The event numbers after selection cuts are given without and with 
beam polarization. Pol+ = Pe-/Pe+ for Pe~=+80% and Pe+=-60%. 



4.3 Contributions to the Mass Uncertainties 



In Table [3] is summarised the contributions to the mass uncertainties 



Error source for Y 


Cut-based analysis 


Iterative Discriminant Analysis 


Detector effects(systematics) 
Charm fragmentation (systematics) 
Stop fragmentation(systematics) 
Neutralino Mass(systematics) 
Background Contribution{systematics) 
Sum of experimental systematics 


0.9% 
0.6% 
0.7% 
0.8% 
0.8% 

1.7%(Ami, =0.10 GeV) 


2.4% 
0.5% 
0.7% 
2.2% 
0.1% 

3.4%(Amj,=0.21 GeV) 


Statistical 

Sum of experimental errors 


3.1%(Amt^=0.19GeV) 
3.5%(Ami^=0.24 GeV) 


2.7%(Amtj=0.17GeV) 
4.3%(Amtj=0.28GeV) 


Theory for signal cross-section 


5.5% 


5.5% 


Total error AY 


6.5% {Ami^ =0.42 GeV) 


7.0%(Amj, =0.44 GeV) 



Table 3: Combination of statistical and systematic errors for the determination of the stop 
mass from a threshold-continuum cross-section measurement. In parenthesis is given the 
overall error on the measured mass. An beam spectrum error Am^^^O.J GcV, is included. 



The next to next to leading order (NNLO) QCD corrections are expected to be of the 
same order than the NLO. This is based on the top quark results. Assuming a factor two 
improvement in the calculations by the time ILC is running (A 1% NNLO correction is also 
included for the EW componant). The relic dark matter density is shown below 

5 Conclusions 

We deal with more realistic data, including quarks hadronization and fragmentation, and 
with a lower integrated luminosity, almost by a factor two, but still manage to improve the 
stop mass precision by a factor three comparatively to The results of the dark Matter 



LCWS/ILC2007-June2,C.Milstene. 



1 — \ — I — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — r 




J \ I \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ L 

121 122 123 124 
% 

Figure 4: Dark matter Relic Density 

relic density versus the stop mass precision are shown in three cases, in the last figure. The 
light grey dots represent our previous results[2] for Amf^= 1.2 GeV, the dark grey dots 
correspond Amf^= 0.42 GeV, QcDAth^ = 0.109+0.015-0.013, include both experimental 
and theoretical errors. The black dot imply 6171^^= 0.24 GeV . ncDMh"^ = 0.109+0.0012- 
0.0010, experimental errors sequential cuts. It is only a small improvement in the precision 
of the dark matter density with respect to Amj^=0.42. The red star is our working point. 
The precision is very comparable to 0.103< flcDMh^ <0.116, the current WMAP results 
are shown by the horizontal green bands on the figure for la and 2a constraints. 
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